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EFFECTING SAME, filed October 15, 1998. 



BACKGROUND OF THE INVENTION 
Technical Field : This invention relates to the field of signal processing. More 

15 particularly, this invention relates to processing measured signals to remove unwanted 
signal components caused by noise and especially noise caused by motion artifacts. 

State of the Art : The measurement of physiological signals can often be difficult 
because the underlying physiological processes may generate very low level signals. 
Furthermore, interfering noise is inherent in the body and the interface between the body 

20 and sensors of physiological processes. Examples of physiological measurements 
include: measurement of electrocardiogram (ECG) signals based on the electrical 
depolarization of the heart muscle, blood pressure, blood oxygen saturation, partial 
pressure of C0 2 , heart rate, respiration rate, and depth of anesthesia. ECG signals, for 
example, are typically detected by surface electrodes mounted on the chest of a patient. 

25 ECG signals are weak at the signal source (i.e., the heart) and are even weaker at the 
surface of the chest. Furthermore, electrical interference from the activity of other 
muscles (e.g., noise caused by patient breathing, general movement, etc.) causes 
additional interference with physiological signals such as an ECG. Thus, considerable 
care must be taken in the design and use of physiological processors to enhance the 

30 quality of the true signal and reduce the effects of interfering noise signals. 
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It is convenient to characterize a measured signal as being a composite signal 
composed of a true signal component and a noise signal component. The terms 
"measured signal" and "composite signal" will be used interchangeably hereinafter. 
Signal processors are frequently used to remove noise signal components from a 
5 composite measured signal in order to obtain a signal that closely, if not identically, 

represents the true signal. Conventional filtering techniques such as low pass, band pass, 
and high pass filtering can be used to remove noise signal components from the measured 
composite signal where the noise signal component occupies a frequency range outside 
the true signal component. More sophisticated techniques for conventional noise filtering 
10 include multiple notch filters, which are suitable for use where the noise signal 

component exists at multiple, distinct frequencies, all outside the true signal frequency 
band. 

However, it is often the case that the frequency spectrum of the true and noise 
signal components overlap and that the statistical properties of both signal components 

15 change with time. More importantly, there are many cases where little is known about 
the noise signal component. In such cases, conventional filtering techniques may be 
ineffective in extracting the true signal. 

The measurement of oxygen saturation in the blood of a patient is a common 
physiological measurement, the accuracy of which may be compromised by the presence 

20 of noise. Knowledge of blood oxygen saturation can be critical during surgery. There are 
means of obtaining blood oxygen saturation by invasive techniques, such as extracting 
and testing blood removed from a patient using a co-oximeter. But, such invasive means 
are typically time consuming, expensive, and uncomfortable for the patient. Fortunately, 
non-invasive measurements of blood oxygen saturation can be made using known 

25 properties of energy attenuation as a selected form of energy passes through a bodily 
medium. Such non-invasive measurements are performed routinely with a pulse 
oximeter. 

The basic idea behind energy attenuation measurements as employed in pulse 
oximetry is as follows. Radiant energy is directed toward a bodily medium, where the 
30 medium is derived from or contained within a patient, and the amplitude of the energy 
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transmitted through or reflected from the medium is then measured. The amount of 
attenuation of the incident energy caused by the medium is strongly dependent on the 
thickness and composition of the medium through which the energy must pass, as well as 
the specific form of energy selected. Information about a physiological system can be 
5 derived from data taken from the attenuated signal of the incident energy transmitted or 
reflected. However, the accuracy of such information is reduced where the measured 
signal includes noise. Furthermore, non-invasive measurements often do not afford the 
opportunity to selectively observe the interference causing the noise signal component, 
making it difficult to remove. 

10 A pulse oximeter is one example of a physiological monitoring system that is 

based upon the measurement of energy attenuated by biological tissues and substances. 
More specifically, a pulse oximeter measures the variable absorption caused by blood 
volume changes, primarily arterial in origin. Pulse oximeters transmit electromagnetic 
energy at two different wavelengths, for example at 660 nm (red) and 940 nm (infrared, 

1 5 hereinafter IR) into the tissue and measure the attenuation of the energy as a function of 
time. The output signal of a pulse oximeter is sensitive to the pulsatile portion of the 
arterial blood flow and contains a component that is a waveform representative of the 
patient's arterial pulse. This type of signal, which contains a component related to the 
patient's pulse, is called a plethysmographic waveform or plethysmogram. 

20 The period of rhythmic contraction of the heart by which blood is driven through 

the aorta and pulmonary artery is known as systole. Maximum light absorbance occurs 
during the systole of a cardiac cycle and is indicated on a plethysmogram by a low point 
or systolic valley. Conversely, the period of rhythmic relaxation and dilation of the heart 
cavities occurs during diastole when blood is drawn into the heart cavities. Minimum 

25 light absorbance occurs during the diastole of a cardiac cycle and is indicated on a 
plethysmogram by a high point or diastolic peak. 

Pulse oximetry measurements typically use a digit, such as a finger, or an ear lobe 
or other element of the body, where blood flows close to the skin as the medium through 
which light energy is transmitted. The finger, for example, is composed of various 

30 tissues and substances including skin, fat, bone, muscle, blood, etc. The extent to which 
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each of these biological tissues and substances attenuate incident electromagnetic energy 
is generally known. However, the effect of motion can cause changes in the optical 
coupling of the sensor (or probe) to the finger, the underlying physiology, the local 
vasculature, optical properties of tissues due to changing optical path length as well as 
5 combinations and interactions of the all of the above. Thus, patient motion may cause 
erratic energy attenuation. 

A typical pulse oximeter includes a sensor, cabling from the sensor to a computer 
for signal processing and visual display, the computer and visual display typically being 
included in a patient monitor. The sensor typically includes two light emitting diodes 
10 (LEDs) placed across a finger tip and a photodetector on the side opposite the LEDs. The 
detector measures both transmitted light signals once they have passed through the finger. 
The signals are routed to a computer for analysis and display of the various parameters 
measured. 

The underlying physical basis of a pulse oximeter is Beer's law (also referred to 
15 as Beer-Lambert's or Bouguer's law) that describes attenuation of monochromatic light 
traveling through a uniform medium that absorbs light with the equation: 

transmitted incident ' \ X J 

20 where Itransmitted is the intensity of the light transmitted through the uniform medium, 
incident is the intensity of incident light, d is the distance light is transmitted through the 
uniform medium, c is the concentration of the absorbing substance in the uniform 
medium, expressed in units of mmol L" 1 , and a(X) is the extinction or absorption 
coefficient of the absorbing substance at wavelength X, expressed in units of L/(mmol 

25 cm). The properties of Beer's law are valid even if more than one substance absorbs light 
in the medium. Each light absorbing substance contributes its part to the total 
absorbance. However, Beer's law does not strictly apply since an LED's output is not 
monochromatic and scattering effects do have a significant influence. Thus, 
manufacturers often utilize an empirically determined lookup table to map from the ratio 

30 of absorbance (or transmittance) at the red and IR frequencies to a saturation value. 
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Two LEDs emit narrowband light (i.e. 9 half power bandwidth of typically 15 nm) 
at two different frequency bands, typically red (centered at about 660 nm) and IR 
(centered at about 940 nm). The intensity of light transmitted through tissue, 
^transmitted* is different for each wavelength of light emitted by the LEDs. 
5 Oxyhemoglobin (oxygenated blood) tends to absorb IR light, whereas deoxyhemoglobin 
(deoxygenated blood) tends to absorb red light. Thus, the absorption of IR light relative 
to the red light increases with oxyhemoglobin. The ratio of the absorption coefficients 
can be used to determine the oxygen saturation of the blood. 

To estimate blood oxygen saturation, SpC>2, a two-solute concentration is 
10 assumed. A measure of functional blood oxygen saturation level, Sp02, can be defined 
as: 



Sp0 2 = 100 — (2) 
c,+c 0 

15 where cq represents oxyhemoglobin solute concentration, and c r represents reduced or 
deoxyhemoglobin solute concentration. 

Noise signal components in a measured pulse oximetry light signal can originate 
from both AC and DC sources. DC noise signal components may be caused by 
transmission of electromagnetic energy through tissues of relatively constant thickness 

20 within the body, e.g., bone, muscle, skin, blood, etc. Such DC noise signal components 
may be easily removed with conventional filtering techniques. AC noise signal 
components may occur when tissues being measured are perturbed and, thus, change in 
thickness while a measurement is being made. Such AC noise signal components are 
difficult to remove with conventional filtering techniques. Since most materials in and 

25 derived from the body are easily compressed, the thickness of such matter changes if the 
patient moves during a non-invasive physiological measurement. Thus, patient 
movement can cause the properties of energy attenuation to vary erratically. The erratic 
or unpredictable nature of motion artifacts induced by noise signal components is a major 
obstacle in removing them. 
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Various approaches to removing motion artifacts from measured physiological 
signals, and particularly for use in pulse oximeters, have been proposed. U.S. Patent Nos. 
5,482,036, 5,490,505, 5,632,272, 5,685,299, 5,769,785 and 6,036,642, all to Diab et al, 
and U.S. Patent No. 5,919,134 to Diab, disclose methods and apparatuses for removing 
5 motion artifacts using adaptive noise cancellation techniques. The basic proposition 
behind these Diab et al. patents is to first generate a noise reference signal from the two 
measured signals, and then use the noise reference signal as an input to an adaptive noise 
canceller along with either or both of the measured signals to remove the reference noise 
signal from the measured signals, thus approximating the actual parametric signals of 

1 0 interest. These Diab et al. patents appear to require the use of both measured input 
signals to generate a noise reference signal. Where the adaptive noise cancellation 
involves the use of a correlation canceller as disclosed in U.S. Patent No, 5,482,036, 
addition problems include significant computational overhead and under certain 
circumstances, the correlation canceller will drive the output signal to zero. 

1 5 Another approach to noise artifact elimination is disclosed in U.S. Patent No. 

5,588,427 to Tien. Tien uses fractal dimension analysis to determine the complexity of 
waveforms in order to determine the proper value of the ratio of true intensities based on 
signal complexity. The Tien approach employs a fractal analyzer to determine values for 
two ratios, a and (3, based on the measured time varying intensity of the transmitted red 

20 and IR light signals including noise, a is defined as the ratio of the time varying true 
intensity of light transmitted from the red LED and the time varying true intensity of the 
light transmitted from the IR LED. p is a similar ratio relating the noise introduced 
during the measurement of the light transmitted by the red LED and the noise introduced 
during the measurement of the light transmitted by the IR LED. According to Tien, a 

25 fractal analyzer then determines values for a and p and provides (a,p) pairs to a statistical 
analyzer. The statistical analyzer performs analysis of one or more (ce,p) pairs to 
determine the best value for a, which is then provided to a look-up table. The look-up 
table provides a value corresponding to the arterial oxygen saturation in the patient. 
While the Tien approach appears to be an innovative use of fractal analysis, it also 

30 appears to be computationally complex. 
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Yet another approach to noise artifact elimination is disclosed in U.S. Patent Nos. 
5,885,213, 5,713,355, 5,555,882 and 5,368,224, all to Richardson et al. The basic 
proposition behind the Richardson et al. approach is to switch operative frequencies 
periodically based on evaluating the noise level associated with various possible 
5 frequencies of operation in order to select the frequency of operation that has the lowest 
associated noise level. It would appear that data measured at a noisy frequency, using the 
Richardson et al. approach could be invalid or useless for calculating arterial oxygen 
saturation. Furthermore, Richardson et al. requires a computational overhead to 
constantly monitor which frequency of operation provides the least noise. 

1 0 Another approach to noise artifact elimination is disclosed in U.S. Patent No. 

5,853,364 to Baker, Jr. et al. The Baker, Jr. et al. approach first calculates the heart rate 
of the patient using an adaptive comb filter, power spectrum and pattern matching. Once 
the heart rate is determined, the oximetry data is adaptively comb filtered so that only 
energy at integer multiples of the heart rate are processed. The comb filtered data and the 

1 5 raw oximetry data are filtered using a Kalman filter to adaptively modify averaging 

weights and averaging times to attenuate motion artifact noise. The adaptive filtering of 
the Baker, Jr. et al. approach appears to add significant computational complexity to solve 
the problem of motion artifact rejection. 

Still another approach to noise artifact elimination is disclosed in U.S. Patent No. 

20 5,43 1 , 1 70 to Mathews. Mathews couples a conventional pulse oximeter light transmitter 
and receiver with a transducer responsive to movement or vibration of the body. The 
transducer provides an electrical signal varying according to the body movements or 
vibrations, which is relatively independent of the blood or other fluid flow pulsations. 
Mathews then provides means for comparing the light signals measured with the 

25 transducer output and performing adaptive noise cancellation. An apparent disadvantage 
of the Mathews approach is the need for a secondary sensor to detect motion. 

Still yet another approach to noise artifact elimination is disclosed in U.S. Patent 
No. 6,002,952 to Diab et al (hereinafter the '952 patent). Diab et al. recognizes the 
limitations of adaptive noise cancellation and particularly the use of a correlation 

30 canceller. The 6 952 patent discloses the use of frequency domain analysis to extract a 
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pulse rate from oximetry data. According to the '952 patent, coupling coefficients related 
to ratios of uncontaminated measurement data and contaminated (noisy) measurement 
data can be determined from taking the ratios at each of a series of spectral peaks 
identified in the frequency domain. The c 952 patent further discloses using the coupling 
5 coefficients to identify the presence of noise by calculating the difference between the 
largest and smallest ratio lines for all spectral peaks, determining whether that difference 
is greater than a pre-selected threshold and whether the frequencies associated with the 
largest and smallest spectral peaks are arbitrarily close or not to each other. Where noise 
is detected, a scale factor is used to scrub the measurement data by controlling the gain 

1 0 control input of a gain controlled amplifier. The scale factor is zero in the presence of no 
noise, and can range up to the largest ratio line where there is noise and the frequencies 
are not close together. However, the signal scrubbing disclosed in the '952 patent 
appears to rely on a very limited measure of noise, i.e., whether the difference between 
the largest and smallest ratio lines is greater than a pre-selected threshold and how close 

1 5 the associated frequencies of largest spectral peak and the smallest spectral peak are 

relative to one another. It would be preferable to have multiple confidence measures in a 
method or system for determining physiological parameters in the presence of motion 
artifacts, e.g., a robust pulse oximeter. 

Thus, a need in the art exists for a method, apparatus and system to eliminate 

20 motion-induced noise artifacts from light signals, that is relatively simple 

computationally, and that does not require more than one sensor, does not use correlation 
cancellers or adaptive noise cancellation and that uses multiple measures of confidence to 
determine physiological parameters accurately. 

25 SUMMARY OF THE INVENTION 

The present invention includes methods, apparatuses and systems for removing 
noise in physiological measurements caused by motion or other similar artifacts. The 
methods, apparatuses and systems of the present invention eliminate noise from light 
signals using a single conventional sensor and are relatively simple computationally. 

30 In accordance with one aspect of the invention, a method of removing motion 
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artifacts from electrical signals representative of attenuated light signals, includes 
transforming the electrical signals into frequency domain data, identifying a plurality of 
candidate peaks from the frequency domain data, analyzing each of the plurality of 
candidate peaks in the context of selected parameters calculated with respect thereto and 
5 arbitrating between each of the plurality of candidate peaks based on the selected 
parameters to select a best frequency. 

In accordance with another aspect of the invention, a method of determining pulse 
rate and saturation from electrical signals representative of attenuated light signals and 
motion artifacts, includes acquiring a segment of red data and a segment of IR data from 

10 each of the electrical signals representative of attenuated light signals, transforming both 
the segment of red data and the segment of IR data into red and IR frequency domain 
data, respectively, identifying a plurality of candidate peaks from the red and IR 
frequency domain data, analyzing each of the plurality of candidate peaks in the context 
of selected parameters calculated with respect thereto, arbitrating between each of the 

15 plurality of candidate peaks based on the selected parameters to select a best frequency, if 
one exists, outputting pulse rate and saturation from the best frequency, and repeating the 
above steps for new segments of data. Additionally, various quality or confidence 
measures may be used to evaluate the validity of the candidates. 

A circuit card embodiment includes a processor with memory for storing a 

20 computer program that is capable of executing instructions embodying methods of the 
invention. 

A system embodiment includes an input device, an output device, a memory 
device and a motion artifact rejection circuit card capable of executing instructions stored 
in the memory device implementing the methods described herein. 
25 Finally, a system embodiment includes an input device, and output device, a 

memory device and a processor, which may be a digital signal processor, capable of 
executing instructions stored in the memory device implementing the methods described 
herein. 
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BRIEF DESCRIPTION OF THE DRAWING FIGURES 
In the drawings, which illustrate what is currently regarded as the best mode for 
carrying out the invention and in which like reference numerals refer to like parts in 
different views or embodiments: 
5 FIG. 1 is a high-level flowchart of a method embodiment of the invention. 

FIG. 2 is two of graphs showing acquired IR and red data segments. 
FIG. 3 is a graph of the power spectrum of the IR data segment in FIG. 2 in 
accordance with the invention. 

FIG. 4 illustrates example graphs of measured IR and red data segments in 
1 0 accordance with the invention. 

FIG. 5 is a graph of the frequency domain transformed IR signal from FIG. 4. 
FIG. 6 illustrates three graphs of IR data after filtering with three different IR 
filters and segmented with vertical lines to show pulses and parameter calculations 
according to the invention. 
15 FIG. 7 is a block diagram of a motion artifact rejection circuit card configured to 

remove noise artifacts from signals representing bodily parameters in accordance with the 
invention. 

FIG. 8 is a block diagram of a pulse oximetry system including a motion artifact 
rejection circuit card capable of removing noise from pulse oximetry data in accordance 
20 with the invention. 

FIG. 9 is a block diagram of a pulse oximetry system including a processor device 
programmed to remove noise from pulse oximetry data in accordance with the invention. 

DETAILED DESCRIPTION OF THE INVENTION. 
25 This patent application is a continuation-in-part of U.S. patent application serial 

number 09/410,991, titled METHOD, APPARATUS AND SYSTEM FOR REMOVING 
MOTION ARTIFACTS FROM MEASUREMENTS OF BODILY PARAMETERS, filed 
October 1, 1999, the contents of which are expressly incorporated herein by reference. 

The following detailed description discloses methods, apparatuses and systems for 
30 removing motion artifacts from measured plethysmography waveforms, particularly, but 
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without limitation, those used in pulse oximetry. A system embodiment of the invention 
includes pulse oximetry hardware and associated software to perform the motion artifact 
suppression. A method embodiment of this invention includes a series of steps that 
exploit certain characteristics of plethysmographic waveforms. The methods, apparatuses 
5 and systems described below are suitable for use with sensors employing light 

transmitted or reflected through bodily tissues and substances. For convenience, the 
following detailed description will assume measurement of light that has been transmitted 
through a finger of a human. The terms "signal" and "waveform" are used 
interchangeably herein. 

10 FIG. 1 is a high-level flowchart of an embodiment of a method of removing 

motion artifacts from plethysmographic data and obtaining a measure of pulse rate and 
Sp0 2 from that data. The method steps include acquiring segments of raw 
plethysmographic data 100, both a red data segment and an IR data segment, conditioning 
each segment of raw data for signal processing 110, transforming the conditioned data 

15 into the frequency domain 120, analyzing the frequency domain data for candidate 

spectral peaks 130, calculating selected parameters associated with the candidate spectral 
peaks 140, arbitrating between the candidate peaks based on the selected parameters to 
select a best frequency 150, outputting pulse rate and median Sp0 2 for the best frequency, 
if a best frequency was found 160, and repeating these steps for new raw data segments 

20 170, as required. The method embodiment of the invention is applied to both red and IR 
data signals to eliminate or reduce noise from the data signals prior to outputting pulse 
rate and Sp0 2 . In the preferred embodiment of the invention, both pulse rate and median 
Sp0 2 are output for valid best frequencies. 

The method of this invention begins with acquiring a segment of data (e.g., five or 

25 more pulses or approximately ten seconds) measured from a single light source 

transmitted through a finger and detected with a sensor on the opposite side of the finger. 
Acquiring a data segment is depicted by block 100 of FIG. 1. FIG. 2 illustrates sample 
segments of IR and red data acquired according to block 100 of FIG. 1. The horizontal 
axis of FIG. 2 is measured in units of time, and specifically here in seconds. The vertical 

30 axis of FIG. 2 is measured in arbitrary units, and specifically here in analog-to-digital 
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output units. For convenience, a 10.24 second segment of data will be used to illustrate 
the method. A 10.24 second segment of data corresponds to 1024 data points with a 
sampling rate of 100 data points per second. It should be readily apparent to one of 
ordinary skill in the art that the method of the invention is not limited to data segments of 
5 this size. The signal processing steps described herein may be performed on both red and 
IR data segments independently and simultaneously. Thus, while the steps of the method 
may be illustrated with data from an IR light signal, the same steps are applicable to data 
from a red light signal and vice versa. The terms "data segment", "input waveform", 
"data signal" and "signal" are used interchangeably herein. 

1 0 A segment of data may be received from a sensor that converts transmitted or 

reflected light signals into electrical signals. U.S. Patent Nos. 5,190,038, 5,398,680, 
5,448,991 and 5,820,550 to Poison et al, the disclosures of each of which are 
incorporated herein by reference, disclose and claim electronic systems for receiving red 
and IR data from a sensor, pre-conditioning the electrical signals and then converting the 

1 5 pre-conditioned electrical signals into digital data using an analog-to-digital converter for 
subsequent digital signal processing. The raw red and IR waveforms may be sampled at 
any convenient data rate. However, for simplicity of illustration, a sampling rate of 100 
Hz will be assumed. Additionally, pulse rate and Sp0 2 may be calculated on any 
convenient periodic or non-periodic basis. However, again for simplicity, we will assume 

20 that pulse rate and Sp0 2 are calculated on a periodic basis every V% second. 

Once a segment of data from a single electrical signal (i.e., Red or IR) has been 
acquired and digitized, it may be conditioned for subsequent signal processing as 
depicted by block 110 of FIG. 1. Signal conditioning may include filtering to reduce 
spectral leakage resulting from subsequent frequency analysis. There are several window 

25 filters that may be suitable for such purposes. For example, and not by way of limitation, 
a Hanning window may be used to reduce spectral leakage. It will be readily apparent to 
one of ordinary skill in the art of digital signal processing that other window filters and 
methods of filtering data to reduce spectral leakage may be selected. As such methods of 
filtering and various filters are known to one of ordinary skill in the art of signal 

30 processing, they will not be further detailed herein. FIG. 3 illustrates the power spectrum 
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of the IR data segment of FIG. 2 after filtering. The vertical axis of FIG. 3 may be 
measured in any arbitrary units of power. The horizontal axis is measured in any units of 
frequency, specifically here in units of bpm. 

The conditioned data is then transformed into the frequency domain for further 
5 analysis and signal processing, see block 120 of FIG. 1. Signal processing as described 
herein is generally performed in the frequency domain. The segment of data is converted 
into the frequency domain by, for example, performing the conventional Fast Fourier 
Transform (FFT) on the data segment. FIG. 3 is a graph of the FFT of the IR data 
segment of FIG. 2. FIG. 3 illustrates a primary candidate peak at a frequency of 

1 0 approximately 70 bpm and a secondary candidate peak at a frequency of approximately 
128 bpm. Other common techniques of converting time-domain data to the frequency 
domain may also be used, e.g., classical methods using the FFT such as the periodogram 
or correlogram, autoregressive methods, Prony's method, minimum variance methods, 
maximum likelihood methods. Additionally, time domain data may be converted to the 

15 frequency domain using transforms such as discrete cosine transform, wavelet transform, 
discrete Hartley transform, and Gabor transform. The preferred transform according to 
this method is the FFT with a window size of 1024 points. The 1024 data points are 
placed in a buffer, the FFT buffer. The FFT transforms the 1024 points of data from the 
time domain into the frequency domain. The output of the FFT is 512 points of real and 

20 512 points of imaginary data in the frequency domain. From these 5 1 2 points of real and 
512 points of imaginary data the power spectrum is calculated and placed in a power 
spectrum buffer. 

Both transient and periodic noise artifacts can induce peaks in the frequency 
domain that may be larger than the peak associated with the patient's heart rate. The 
25 frequency peak that actually represents the patient's heart rate (best frequency) must then 
be determined. Analyzing the power spectrum peaks to determine candidate spectral 
peaks is depicted in block 130 of FIG. 1 . One approach to determining the best frequency 
would be to order the frequencies by peak amplitude from largest to smallest, F, to F n , 
where ¥ x through F n are not harmonics of each other, and analyze them one by one to find 
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the correct frequency, i.e., the patient's heart rate. However, the preferred method selects 
up to three candidate spectral peaks for further analysis. 

The function of block 130 is to locate candidate spectral peaks from the power 
spectrum computed in block 120. The power spectrum buffer is an array of 512 vector 
5 points (referred to herein as "bins") in the frequency domain. Each array element in the 
power spectrum buffer represents the power of the corresponding frequency in the 
original raw data waveform. Of the 512 bins, only bins 5 (29 bpm) through 43 (252 bpm) 
are of interest, since this range covers the physiological limits of the human heart rate. 
All other bins are unused by the method of the invention because they cannot 
10 physiologically represent a valid spectral frequency of a pulse rate. Table 1, below, 
shows the first 45 points of the power spectrum array. 



Table 1 



Power Spectrum Buffer 


Frequency (Hz) 


Pulse Rate (bpm) 


bin number n 


f=n* 100/1024 


Pulse Rate =/* 60 


0 


0.00000 


0.0 


1 


0.09766 


5.9 


2 


0.19531 


11.7 


3 


0.29297 


17.6 


4 


0.39063 


23.4 


5 


0.48828 


29.3 


6 


0.58594 


35.2 


7 


0.68359 


41.0 


8 


0.78125 


46.9 


9 


0.87891 


52.7 


10 


0.97656 


58.6 


11 


1.07422 


64.5 


12 


1.17188 


70.3 


13 


1.26953 


76.2 


14 


1.36719 


82.0 


15 


1.46484 


87.9 


16 


1.56250 


93.8 


17 


1.66016 


99.6 


18 


1.75781 


105.5 


19 


1.85547 


111.3 
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20 


1.95313 


117.2 


21 


2.05078 


123.0 


22 


2.14844 


128.9 


23 


2.24609 


134.8 


24 


2.34375 


140.6 


25 


2.44141 


146.5 


26 


2.53906 


152.3 


27 


2.63672 


158.2 


28 


2.73438 


164.1 


29 


2.83203 


169.9 


30 


2.92969 


175.8 


31 


3.02734 


181.6 


32 


3.12500 


187.5 


33 


3.22266 


193.4 


34 


3.32031 


199.2 


35 


3.41797 


205.1 


36 


3.51563 


210.9 


37 


3.61328 


216.8 


38 


3.71094 


222.7 


39 


3.80859 


228.5 


40 


3.90625 


234.4 


41 


4.00391 


240.2 


42 


4.10156 


246.1 


43 


4.19922 


252.0 


44 


4.29688 


257.8 



In table 1, column 1 is the bin number, n; column 2 is the center frequency,/ of 
the corresponding bin number, w, calculated as the product of the bin number and 
sampling rate (100 samples/sec) divided by the block size used by the FFT (i.e., 1024); 
5 and column 3 is the pulse rate corresponding to the center frequency,/ of column 2, 
calculated by multiplying/(measured in units of beats per second) by 60 to convert to 
units of beats per minute. 

In order to select candidate peaks (and corresponding frequencies), different 
amplitude analysis methods are applied to different frequency bands. The amplitude of 
10 adjacent and nearby frequency components of the candidate peak amplitude may be 
compared in terms of their absolute or relative values. For example, the frequencies 
represented by candidate bins 5 through 10 ("5-10" search method) may be stepped 
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through in a sequential fashion. According to the "5-10" search method, a candidate bin 
is assumed to be a candidate power spectrum peak if the previous three bins and 
subsequent four bins relative to a candidate bin are all lower in power than the candidate 
bin. For example, in order for bin 6 to be a candidate spectral peak, bins 3, 4, 5, 7, 8, 9 
5 and 10 must all be lower in power than bin 6. The terms "spectral peak", "power peak", 
or simply "peak" are used synonymously herein. Various amplitude, shape, syntactic or 
other pattern analysis methods may be applied to identify a candidate peak. Also, 
multiple curve fit methods, as known to one of ordinary skill in the spectroscopic 
analysis, may also be applied. 

10 Once all possible power spectrum peak candidates are found, predetermined 

criteria are applied to select, at most, three candidate spectral peaks. First, the power 
peak associated with the largest power amplitude is selected to be the primary candidate 
peak. Then, any power peaks that are determined to be harmonics of the primary 
candidate power peak are eliminated. According to the method, a harmonic is defined as 

1 5 any power peak the frequency of which is a multiple of the primary peak, ± 1 bin, and the 
amplitude of which is less than half the maximum allowed power of the previous 
harmonic, or in the case of the first harmonic, less than half the power of the primary 
peak. For example, assume a primary candidate peak is found at bin 10. Possible 
harmonic bins of 10 are bins 19-21,29-31 and 39-41. Continuing with the example, if 

20 the primary power peak amplitude (bin 1 0) is 1 00 arbitrary power spectrum units, then 
bins 19-21 must be less than 50 units to be deemed a harmonic, bins 29-31 must be less 
than 25 units and bins 39-41 must be less than about 12 units, where units are the 
measure of the amplitude of the power spectrum. Other weights may be applied to the 
analysis of the sequence for detection of harmonics of the candidate spectral peak without 

25 departing from the scope of the invention. 

After harmonics of the primary candidate peak are eliminated, the next largest 
remaining power peak found (if any) is selected to be the secondary candidate peak. 
Finally, if the previous pulse rate is non-zero, the power spectrum corresponding to the 
previous pulse rate is determined. If the bin corresponding to the previous pulse rate is 

30 not equal to the primary or secondary candidate power peak, then the bin corresponding 
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to the previous pulse rate is selected to be the tertiary candidate peak. Thus, up to three 
candidate peaks (primary, secondary and tertiary) and corresponding frequencies of each 
candidate peak are identified in block 130 of FIG. 1. 

Block 140 of FIG. 1 depicts calculating selected parameters associated with the 
5 candidate peaks identified in block 130. Block 140 may include pulse window filtering, 
and calculating such parameters as peak detection, pulse rejection criteria and descriptive 
parameters associated with each of the up to three candidate power peaks found from 
block 130. These parameters are used to determine a pulse confidence for each candidate 
peak. The parameters calculated according to the invention for each filtered candidate 

10 peak include measures of central tendency and variability of pulse width, pulse rate and 
Sp0 2 , as well as measures of the history and confidence of these parameters. The 
preferred embodiment includes parameters such as: (1) Window Pulse Rate, (2) Pulse 
Width Variability, (3) Sp02 Variability, (4) Pulse Window Sp02, (5) Pulse Peak 
Amplitude Variability, (6) Pulse Rate History Percentage, and (7) Pulse Window 

15 Confidence. It should be noted that other parameters including, but not limited to, other 
measures of central tendency, variability (i.e., skewness, kurtosis), history/trend and 
confidence, could be calculated from the candidate power peaks without departing from 
the scope of the invention. Each of the parameters listed is discussed in greater detail 
below, beginning with pulse window filtering. 

20 Prior to calculating the aforementioned parameters, each candidate peak is filtered 

with a narrow band filter, such as a bandpass, finite impulse response (FIR) filter. In one 
aspect of the present invention, one of several predefined FIR filters is applied to a given 
bin or candidate peak. The peak frequencies of the filters may be separated by a fixed 
difference in frequency (measured in Hz or bpm), such as 25 bpm, or may be variable and 

25 a function of either frequency or a characteristic of the spectrum, for example variability 
or noise, or both. For example, if a candidate peak was found at bin 12, a filter with 
center or peak frequency of 76.2 bpm might be chosen. A fixed difference in frequency 
may be in a range from about 15 bpm to about 40 bpm. Likewise, a variable difference in 
frequency may be in a range from about 15 bpm to about 40 bpm. 
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Preferably, to improve discrimination, especially with closely spaced peaks, the 
bandpass filter coefficients may be stored or generated and adjusted as needed so that the 
center frequency is nearly identical to the candidate frequency. Additionally, other 
filtering methods such as (a) other types of bandpass filters, i.e., infinite impulse response 
5 (IIR) filters, and (b) frequency domain methods such as transforming the data into the 
frequency domain (for example, FFT), filtering or eliminating unwanted components in 
the frequency domain and transforming back into the time domain (for example, inverse 
FFT) for further signal processing, may be applied. 

Once the up to three peak candidates are selected and filtered, a peak detector 

10 algorithm is applied to each of the up to three candidate peaks in the time domain. The 
function of the peak detector algorithm is to identify power spectrum peaks in each of the 
filtered time domain data segments and their associated center frequencies. The terms 
"window" and "pulse window" are used interchangeably with "time domain data 
segment" herein. For each peak found in the time domain, the pulse width is calculated 

15 as the time between each peak. The Window Pulse Rate is calculated by dividing the sum 
of the pulse width time of all peaks by the number of peaks detected. 

Pulse Width Variability, a measure of how consistent the pulse width is for all the 
peaks in a given pulse window, is calculated according to the method of the invention. 
With the exception of subjects presenting cardiac arrhythmias, particularly ventricular 

20 arrhythmias, the variability of the pulse width of all the peaks should be low within a 
short time interval such as a 10.24 second window. Higher pulse width variability often 
is an indication of either (a) cardiac arrhythmias or (b) physiological artifacts such as 
motion. Pulse Width Variability is calculated as the sum of absolute differences between 
individual pulse widths and the average pulse width normalized by the average pulse 

25 width: 

^ | A verage Pulse Width - Pulse Width, | 

Pulse Width Variability = - 1 , (3) 

Average Pulse Width 
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where, i is the number of peaks detected in the window, Pulse Widthj is the pulse width 
for the z'th peak, and Average Pulse Width is the sum of the individual pulse widths 
divided by the number of pulses. For example, a pulse rate of 1 80 bpm has a pulse width 
of 330 ms; and a pulse rate of 60 bpm has a pulse width of 1000 ms. If an average pulse 
5 width difference was found to be 1 00 ms, this would have a much greater effect at 1 80 
bpm than 60 bpm. Thus, dividing by the pulse rate normalizes the pulse width 
variability. Returning to the 180 bpm versus 60 bpm example, dividing by the Average 
Pulse Width causes the Pulse Width Variability to be 3 times greater at 180 bpm than at 
60 bpm. 

1 0 Sp0 2 is calculated for each peak in the pulse window using the ratio, R, which is 

'mapped 5 to Sp0 2 via a lookup table. This ratio, R, is defined as: 



Red AC Component 
Red DC Component 
f IRAC Component ^ 
K IR DC Component j 



15 The ratio, R, in equation 4 is used to index into an empirically derived table to determine 
Sp0 2 . The IR AC Component is chosen at the point of maximum negative slope between 
the peak and valley for each peak of the filtered IR waveform. The Red AC Component is 
the slope of the filtered Red waveform at the time coincident with the above selected IR 
AC Component. The peak and valley points from the IR filtered waveforms are 

20 transposed onto the raw Red and IR waveform. The average between the peak and the 
valley is considered the DC component (analogous to a DC offset for a positively biased 
AC waveform). This DC component is calculated for both the Red and IR waveforms 
(i e. , Red DC Component and IR DC Component). This process is repeated for each of 
the up to three bandpass filtered pulse windows corresponding to the candidate peaks 

25 identified in block 130 of FIG. 1. The term "pulse window" is used herein to represent 
time domain data corresponding to a particular candidate peak that has been bandpass 
filtered. 
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Sp02 Variability, a measure of how consistent the Sp0 2 is for all the peaks in a 
pulse window, is calculated according to the method of the invention. Under typical 
conditions, SpC>2 Variability is low, often within + 2 percent saturation over a short time 
interval such as the 10.24 second pulse window. When the pulse window is filtered by a 
5 frequency that is not related to the pulse rate (e.g., random noise), the Sp02 Variability 
tends to be high. Therefore, Sp02 Variability is a good measure for determining 
confidence in a pulse window. Sp02 Variability is calculated as the sum of the absolute 
difference between the individual Sp0 2 values and the average Sp0 2 for the pulse 
window. 

10 

Sp0 2 Variability = Y\^verage Sp0 2 - SpQ 2 i | , (5) 



where i is the number of peaks detected in the window, and Sp02, i is the saturation 
calculated for the zth peak detected and Average Sp02 is sum of the individual Sp0 2 
1 5 values divided by the number of individual Sp0 2 values. 

The Pulse Window Sp02 is calculated by the method of the invention as the 
median value of all of the Sp0 2 values within the current pulse window. Other methods 
of determining central tendency may be used including by not limited to a weighted mean 
or average. 

20 The Pulse Peak Amplitude Variability, a measure of the consistency of the 

amplitude of the pulse peaks in a pulse window, is calculated according to the method of 
the invention. Pulse Peak Amplitude Variability is calculated as the sum of the 
differences between the individual pulse peak amplitudes and the average pulse peak 
amplitude for the pulse window. 

25 Motion artifacts are usually not purely rhythmic in nature. Therefore, the portion 

of the power spectrum comprising motion artifacts changes dynamically as the spectrum 
of the motion artifacts changes. In contrast, the spectrum of the underlying pulse rate 
varies much less over longer periods of time relative to motion artifact spectrum. 
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Pulse Rate History Percentage, another parameter useful for detecting motion 
artifacts, is calculated by the method of the invention. According to the method, a pulse 
rate is calculated for the primary and secondary candidate peaks and these pulse rate 
calculations are saved in memory. This memory may be any capacity but preferably is 
capable of storing at least between 10-60 seconds of pulse rates (for the primary and 
secondary candidate peaks) and is updated periodically with the newest values 
overwriting the oldest values. For illustration purposes, assume the memory stores 30 
seconds of pulse rates for the primary and secondary peaks. Pulse Rate History 
Percentage is calculated according to the method of the invention as the percentage that 
the pulse rate corresponding to the candidate peak occurred in a given period of time, e.g., 
the most recent 30 seconds. Of course, one of ordinary skill in the art may recognize that 
a Pulse Rate History Percentage may be calculated in other analogous ways. For 
example, there could be a longer history of pulse rates (i.e., more or less than 30 seconds) 
and it could be weighted or filtered in various manners without departing from the scope 
of the invention. 

A pulse window under evaluation may be rejected from further processing and 
flagged as an invalid pulse window if certain criteria are met. A pulse window under 
evaluation may be checked against the following criterion: 

1 . The number of peaks in the pulse window, z, is two or less (/ <2) 

2. The Window Pulse Rate is zero (i.e., no frequency found). 
Additional criteria related to variability and history may include the following: 

1 . The Pulse Rate History Percentage is less than a percentage of a significant 
portion of the pulse. According to the preferred embodiment of the invention, 
a percentage of a significant portion of the pulse would be in the range from 
about 25% to about 30%. 

2. The Sp02 Variability is greater than "normal" variation in Sp0 2 in either 
absolute or relative terms. According to the preferred embodiment of the 
invention, SpC>2 Variability greater than a threshold ranging from about 3% to 
about 5% is greater than "normal" variation in Sp0 2 . 
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3 . The Pulse Width Variability is greater than a threshold representing excessive 
variations. According to the preferred embodiment of the invention, a 
threshold representing excessive variations may fall within the range of about 
200 to about 400 points for a pulse window of 1024 points. 

4. The Window Pulse Rate differs by more than an excessive amount either in 
absolute or relative terms from the center frequency of the candidate power 
spectrum peak. According to the preferred embodiment of the invention, an 
excessive amount is a threshold greater than about 20 bpm to about 30 bpm 5 or 
about 20% to about 35% of the center frequency, whichever maximum 
threshold is smaller. 

If the pulse window under evaluation meets any of these criteria, then the pulse 
window under evaluation is rejected as invalid and flagged as such. Also, the optimal 
thresholds and values for each of the above criteria may be optionally adjusted by 
methods known to one of ordinary skill in the art, including but not limited to, learning or 
search methods. 

A confidence measure, Pulse Window Confidence, is also calculated according to 
block 140 of FIG. 1. According to the preferred method of the invention, Pulse Window 
Confidence is calculated as a weighted sum of the Pulse Width Variability, Sp02 
Variability, Pulse Amplitude Variability and the Pulse Rate History Percentage 
parameters. The lower the value of Pulse Window Confidence measure, the higher the 
confidence that the candidate peak under evaluation is a valid pulse rate. The Pulse 
Window Confidence, which is a point value without units, is only computed for each of 
the up to three remaining candidate peaks and then passed to the arbitrating step, see 
block 150 of FIG. 1. 

The function of the arbitrating step, block 150 of FIG. 1, is to determine which, if 
any, of the up to three candidate peaks should be accepted. The arbitrating step 150, is 
accomplished by evaluating the calculated parameters, including the confidence or 
quality measures, (/.e., Pulse Window Confidence), of the candidate peaks relative to one 
another. Some of the candidate peaks may already have been flagged as an invalid pulse 
window, and thus, are not evaluated further. If none of the up to three candidate peaks is 
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valid, no new pulse rate or new saturation will be displayed according to the method of 
the invention. Alternatively, if none of the up to three candidate peaks are valid, another 
algorithm (other than the method of the invention) may be employed to determine the 
pulse rate and saturation, see for example U.S. Patent Nos. 5,190,038, 5,398,680, 
5 5,448,991 and 5,820,550 to Poison et al. Arbitration is then conducted among the up to 
three remaining candidate peaks in order to determine which, if any candidate peak, 
should be selected as the best frequency. The arbitration is preferably executed in the 
sequence presented below. 

1. If the primary candidate peak frequency, f] 9 is zero, then there is no 
1 0 valid candidate peak (/. e. , no best frequency) . 

2. If the tertiary candidate peak Pulse Window Confidence is less than the Pulse 
Window Confidence for either the primary candidate peak or the secondary candidate 
peak, then the tertiary candidate peak is the best frequency. Recall that the lower the 
Pulse Window Confidence value, the higher the confidence that the candidate peak is the 

15 true pulse rate. 

3. If the primary candidate peak and the secondary candidate peak have both been 
rejected, then there is no valid candidate peak (i.e., no best frequency). 

4. If the primary candidate peak has not been rejected and the secondary 
candidate peak has been rejected, then the primary candidate peak is the best frequency. 

20 5. If the primary candidate peak has been rejected and the secondary candidate 

peak has not been rejected, then the secondary candidate peak is the best frequency. 

6. If the primary candidate peak Pulse Window Confidence is greater than the 
secondary candidate peak Pulse Window Confidence by a specified threshold, tj, and the 
primary candidate peak Pulse Rate History Percentage is greater than another specified 

25 threshold, then the primary candidate peak is the best frequency. Similar criteria 
apply if the secondary candidate peak Pulse Window Confidence is greater than the 
primary candidate peak. 

7. If the secondary candidate peak frequency, is a rough harmonic of the 
primary candidate peak frequency,//, and the Pulse Window Confidence of the primary 

30 candidate peak is not more than a specified number of points greater than the Pulse 
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Window Confidence of the secondary candidate peak, then accept the primary candidate 
peak. Secondary candidate peak frequency, f2, is a rough harmonic of the primary 
candidate peak frequency, fj 9 if the candidate frequency is within a frequency tolerance of 
approximately ± 10 bpm. Again, similar criteria apply if the secondary candidate peak is 
5 a rough harmonic of the primary candidate peak. 

8. If the Pulse Window Confidence of the primary candidate peak is no more than 
a specified number of points greater than the Pulse Window Confidence of the secondary 
candidate peak, then accept the primary candidate peak. Otherwise, accept the secondary 
candidate peak. 

10 Once a candidate peak has been accepted (as the best frequency) according to the 

arbitrating step 150, the pulse rate and Sp0 2 are calculated for the best frequency and 
output, for example, to a display or monitor, as depicted in block 160 of FIG. 1. The 
steps 100-160 may then be repeated for any new segments of data as depicted in decision 
block 170 of FIG. 1. The above sequence is exemplary only, and not intended to be 

15 limiting. Furthermore, one of ordinary skill in the art will recognize that the various 
criteria selected to evaluate pulse shape may be assigned weights to emphasize relative 
importance. 

FIGS. 4-6 illustrate exemplary graphical results from application of the preferred 
method of the invention. FIG. 4 shows graphs of measured IR and red data segments in 

20 accordance with the invention. FIG. 5 is a graph of the frequency domain transformed IR 
signal from FIG. 4, showing a primary candidate peak at approximately 82 bpm, a 
secondary candidate peak at approximately 105 bpm. 

FIG. 6 illustrates three graphs of IR data after filtering with three different FIR 
filters and segmented with vertical lines to delineate pulses and parameter calculations 

25 according to the invention. The parameter calculations displayed to the right in the 
graphs shown in FIG. 6 are exemplary only, and are not necessary for practicing the 
invention. With respect to those parameter calculations displayed, "PRHist" corresponds 
to Pulse Rate History Percentage; "PR" corresponds to Window Pulse Rate; "PWVar" 
corresponds to Pulse Width Variability; "MSat" corresponds to Pulse Window Sp02; 

30 "SVar" corresponds to Sp(?2 Variability; "PConf ' corresponds to Pulse Window 
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Confidence', "Conf OK" corresponds to an accepted candidate peak or the best frequency, 
and "R:xxxx" corresponds to a notation that the candidate peak under evaluation has been 
rejected for the reason "xxxx", i.e., a parameter calculation has concluded with a rejection 
of the candidate peak all as disclosed herein. 
5 The top graph in FIG. 6 represents the pulse window corresponding to the primary 

candidate peak at approximately 82 bpm as shown in FIG. 5. The middle graph in FIG. 6 
represents the pulse window corresponding to the secondary candidate peak at 
approximately 105 bpm as shown in FIG. 5. The bottom graph in FIG. 6 represents the 
pulse window corresponding to the tertiary candidate peak at approximately 158 bpm as 

1 0 shown in FIG. 5. Note that in this instance, the tertiary candidate at frequency of 1 58 
bpm has the lowest Pulse Window Confidence and no rejections based on calculated 
parameters. Note that FIGS. 4-6 are merely exemplary graphs illustrating sample 
calculations based on actual data obtained from typical pulse oximetry measurements 

The methods described above may be integrated into apparatuses and/or systems 

15 for calculating blood oxygen saturation. Referring to FIG. 7, one apparatus embodiment 
of this invention comprises a motion artifact rejection circuit card 10 with an I/O device 
1 1, a processor 12 and memory 14 for storing a computer programmed algorithm for 
motion artifact rejection as described in the above methods. Processor 12 may be a 
digital signal processor. I/O device 1 1 may be any circuitry that allows communication 

20 to and from external circuitry, for example, and not by way of limitation, bus interface 
circuitry. I/O device 1 1 may include a circuit card edge connector for plugging into a 
pulse oximetry monitor system. Memory 14 may be any solid-state electronic memory 
suitable for storing digital data including, for example, computer code and measurement 
data. 

25 Referring to FIG. 8, the motion artifact rejection circuit card 10 of FIG. 7 may be 

incorporated in a complete pulse oximetry system 16 for eliminating motion-induced 
noise artifacts in electrical signals (as described in the method embodiments above) and 
calculating and displaying physiological parameters, either as a discrete circuit card or as 
part of a larger circuit card, such as a motherboard, controlling other functions of the 

30 pulse oximetry system 16. The pulse oximetry system 16 also includes an input device 
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1 8 and an output device 20. Input device 1 8 may be a pulse oximeter sensor with red and 
IR LED light sources and a photodetector to convert transmitted or reflected light into an 
electrical signal. Output device 20 may be a display device such as a cathode ray tube 
device, liquid crystal display, active matrix display or any other suitable display device 
5 known to one of skill in the art. Alternatively, output device 20 may be a printer for 
producing a permanent or written record such as a laser printer, ink jet printer, thermal 
printer, dot matrix printer or any other suitable printer known to one of skill in the art. 
The pulse oximetry system 16 may be any pulse oximeter that uses the principles of 
operation as described above. A particular pulse oximeter for which the circuit card 

10 embodiment as described above is suitable for use is the Novametrix Medical Systems, 
Inc., Model 520A, Pulse Oximeter. 

Referring to FIG. 9, a block diagram of a pulse oximetry system 22 including a 
processor device 12, an input device 18, an output device 20 and a storage device 24, is 
shown. Input device 1 8 may be a pulse oximeter sensor with red and IR LED light 

1 5 sources and a photodetector to convert transmitted or reflected light into an electrical 
signal. Output device 20 may be a display device such as a cathode ray tube device, 
liquid crystal display, active matrix display or any other suitable display device known to 
one of skill in the art. Alternatively, output device 20 may be a printer for producing a 
permanent or written record such as a laser printer, ink jet printer, thermal printer, dot 

20 matrix printer or any other suitable printer known to one of skill in the art. Storage 
device 24 may be a disk drive, or any kind of solid-state electronic memory device 
suitable for storing digital data including, for example, computer code and measurement 
data. 

Although this invention has been described with reference to particular 
25 embodiments, the invention is not limited to these described embodiments. Rather, it 
should be understood that the embodiments described herein are merely exemplary and 
that a person skilled in the art may make many variations and modifications without 
departing from the spirit and scope of the invention. All such variations and 
modifications are intended to be included within the scope of the invention as defined in 
3 0 the appended claims . 
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What is claimed is: 



CLAIMS 



1 . A method of removing motion artifacts from electrical signals 
5 representative of attenuated light signals, comprising: 

transforming the electrical signals into frequency domain data; 
identifying a plurality of candidate peaks from said frequency domain data; 
developing parameters associated with each of said plurality of identified candidate peaks; 
analyzing each of said plurality of candidate peaks with respect to at least some of said 
1 0 developed parameters; and 

arbitrating between at least some of said plurality of candidate peaks employing at least 
some of said developed parameters to select a best frequency. 



2. The method of claim 1, further comprising conditioning said electrical 
15 signals to reduce spectral leakage prior to said transforming step. 

3. The method of claim 2, wherein said conditioning includes filtering said 
electrical signals. 

20 4. The method of claim 3, wherein said filtering is performed with a Hanning 

window. 

5. The method of claim 1 , wherein said transforming the electrical signals 
into frequency domain data is performed with a fast Fourier transform. 

25 

6. The method of claim 1, wherein said transforming the electrical signals 
into frequency domain data is performed with a technique selected from the group 
comprising periodogram, a correlogram, autoregressive methods, Prony's method, 
minimum variance methods, maximum likelihood methods, a discrete cosine transform, a 

30 wavelet transform, a discrete Hartley transform and a Gabor transform. 
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7. The method of claim 1 , wherein said identifying said plurality of candidate 
peaks comprises: 

assigning a largest power amplitude from said frequency domain data as a primary 
candidate peak; 

5 assigning a next largest power amplitude that is not a harmonic of said primary candidate 
peak as a secondary candidate peak; and 
assigning a previous non-zero pulse rate as a tertiary candidate peak if said previous non- 
zero pulse rate is neither said primary candidate peak nor said secondary 
candidate peak. 

10 

8. The method of claim 1 , wherein said identifying said plurality of candidate 
peaks comprises identifying n peaks, by frequency, Fj to F n? in descending order of peak 
amplitude, where F x through F n are not harmonics of each other. 

15 9. The method of claim 1 , further comprising filtering each of said plurality 

of candidate peaks with a narrow band filter before said developing parameters step. 

10. The method of claim 9, wherein said narrow band filter is a finite impulse 
response filter. 

20 

1 1 . The method of claim 9, wherein said narrow band filter is an infinite 
impulse response filter. 

12. The method of claim 1, further comprising filtering each of said plurality 
25 of candidate peaks with one of n narrow band filters before said developing parameters 

step to mask influence of candidate frequencies not under evaluation. 

13. The method of claim 12, wherein n = 8 and wherein each of said 8 narrow 
band filters is separated by a fixed difference in frequency in a range of approximately 25 

30 bpm to approximately 30 bpm. 
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14. The method of claim 12, wherein each of said n narrow band filters is 
separated by a variable difference in frequency in a range of approximately 25 bpm to 
approximately 30 bpm. 

5 

15. The method of claim 1 , further comprising filtering each of said plurality 
of candidate peaks with a narrow band filter of variable center frequency. 

16. The method of claim 1 5 further comprising filtering each of said plurality 
10 of candidate peaks with a narrow band filter wherein filter coefficients are generated and 

adjusted so that the center frequency of said narrow band filter is approximately the 
frequency associated with each of said candidate peaks. 

17. The method of claim 1, further comprising filtering each of said plurality 
15 of candidate peaks using a fast Fourier transform (FFT), narrow band filter and inverse 

FFT. 

1 8. The method of claim 1 , wherein said analyzing each of said plurality of 
candidate peaks with respect to said developed parameters includes calculating a window 

20 pulse rate. 

1 9. The method of claim 1 , wherein said analyzing each of said plurality of 
candidate peaks with respect to said developed parameters includes calculating pulse 
width variability. 

25 

20. The method of claim 1, wherein said analyzing each of said plurality of 
candidate peaks with respect to said developed parameters includes calculating Sp0 2 
variability. 
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21. The method of claim 1, wherein said analyzing each of said plurality of 
candidate peaks with respect to said developed parameters includes calculating pulse 
window Sp0 2 . 

5 22. The method of claim 1, wherein said analyzing each of said plurality of 

candidate peaks with respect to said developed parameters includes calculating pulse rate 
history percentage. 

23. The method of claim 1, wherein said analyzing each of said plurality of 
1 0 candidate peaks with respect to said developed parameters includes calculating pulse 

window confidence. 

24. The method of claim 23, wherein said calculating pulse window 
confidence includes a weighted sum of pulse width variability, Sp0 2 variability and pulse 

1 5 rate history percentage. 

25. The method of claim 1, wherein said analyzing each of said plurality of 
candidate peaks with respect to at least some of said developed parameters includes 
calculating a window pulse rate, pulse width variability, Sp0 2 variability, pulse window 

20 Sp0 2 , pulse rate history percentage and pulse window confidence. 

26. The method of claim 1 , wherein said arbitrating between at least some of 
said plurality of candidate peaks based on said developed parameters includes applying a 
predetermined criteria to select said best frequency. 

25 

27. The method of claim 1 , wherein said plurality of candidate peaks 
comprises up to three candidate peaks, including a primary candidate peak, a secondary 
candidate peak and a tertiary candidate peak. 
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28. The method of claim 27, wherein said arbitrating between each of said up 
to three candidate peaks includes applying the following criteria using at least some of 
said developed parameters to select said best frequency: 

(a) if the primary candidate peak frequency is zero, then there is no valid 
5 candidate peak; 

(b) if the tertiary candidate peak pulse window confidence is less than the pulse window 

confidence for either the primary candidate peak or the secondary candidate peak, 
then the tertiary candidate peak is said best frequency; 

(c) if the primary candidate peak and the secondary candidate peak have both been 
10 rejected, then there is no valid candidate peak; 

(d) if the primary candidate peak has not been rejected and the secondary candidate peak 

has been rejected, then the primary candidate peak is said best frequency; 

(e) if the primary candidate peak has been rejected and the secondary candidate peak has 

not been rejected, then the secondary candidate peak is said best frequency. 
15 (f) if the primary candidate peak pulse window confidence is greater than the secondary 
candidate peak pulse window confidence by a first threshold, t l5 and the primary 
candidate peak pulse rate history percentage is greater than a second threshold, t^ 
then the primary candidate peak is said best frequency; 

(g) if the secondary candidate peak frequency is a rough harmonic of the primary 
20 candidate peak frequency and the pulse window confidence of the primary 

candidate peak not more than a specified number of points greater than the pulse 
window confidence of the secondary candidate peak, then accept the primary 
candidate peak; and 

(h) if the pulse window confidence of the primary candidate peak is no more than a 
25 specified number of points greater than the pulse window confidence of the 

secondary candidate peak, then said primary candidate peak is said best 
frequency, otherwise, said secondary candidate peak is said best frequency. 

29. The method of claim 28, wherein said arbitration is conducted in the 
30 sequence presented in claim 28. 
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30. A method of determining pulse rate and saturation from electrical signals 
representative of attenuated light signals and motion artifacts, comprising: 

acquiring a segment of red data and a segment of IR data from each of the electrical 
5 signals representative of attenuated light signals; 

transforming both said segment of red data and said segment of IR data into red and IR 

frequency domain data, respectively; 
identifying a plurality of candidate peaks from said red and IR frequency domain data; 
developing parameters associated with each of said plurality of identified candidate 
10 peaks; 

analyzing each of said plurality of candidate peaks with respect to at least some of said 

developed parameters; 
arbitrating between at least some of said plurality of candidate peaks employing at least 

some of said selected parameters to select a best frequency; 
15 outputting pulse rate and saturation relating to said best frequency; and 
repeating the above steps for new segments of data. 

3 1 . The method of claim 30, wherein said transforming includes performing a 
fast Fourier transform. 

20 

32. The method of claim 30, wherein said identifying said plurality of 
candidate peaks comprises: 

assigning a largest power amplitude from said frequency domain data as a primary 
candidate peak; 

25 assigning a next largest power amplitude that is not a harmonic of said primary candidate 
peak as a secondary candidate peak; 
assigning a previous non-zero pulse rate as a tertiary candidate peak if said previous non- 
zero pulse rate is neither said primary candidate peak nor said secondary 
candidate peak. 

30 
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33. The method of claim 30, wherein said developed parameters includes 
window pulse rate calculated by dividing sum of all pulse width times of all peaks in a 
data segment by quantity of peaks detected in said data segment. 

5 34. The method of claim 30, wherein said developed parameters includes 

pulse width variability calculated as the sum of absolute differences between individual 
pulse widths and average pulse width normalized by said average pulse width. 

35. The method of claim 30, wherein said developed parameters includes 
10 Sp0 2 variability calculated as a sum of absolute difference between individual Sp0 2 

values and average Sp0 2 for a given pulse window. 

36. The method of claim 30, wherein said developed parameters includes 
pulse window Sp0 2 calculated by taking a measure of central tendency of all individual 

1 5 Sp0 2 calculations in a given pulse window. 

37. The method of claim 30, wherein said developed parameters includes 
pulse peak amplitude variability calculated as a sum of differences between individual 
pulse peak amplitudes and average pulse peak amplitude for a given pulse window. 

20 

38. The method of claim 30, wherein said developed parameters includes 
pulse rate history percentage calculated as a percentage of time that a pulse rate 
corresponding to a candidate peak has occurred in a given period of time. 

25 39. The method of claim 30, wherein said developed parameters includes 

pulse window confidence calculated as a weighted sum of pulse width variability, Sp0 2 
variability, pulse amplitude variability and pulse rate history percentage. 
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40. The method of claim 30, wherein said arbitrating between at least some of 
said plurality of candidate peaks employing at lease some of said developed parameters 
includes applying selection criteria to select said best frequency. 



5 4 1 . A method of removing motion artifacts from a single electrical signal 

representative of an attenuated light signal, comprising: 
acquiring a segment of data from the electrical signal; 
conditioning said segment of data for signal processing; 
transforming said conditioned segment of data into frequency domain data; 
1 0 analyzing said frequency domain data to identify candidate peaks; 
developing parameters associated with each of said candidate peaks; 
arbitrating between said candidate peaks using at least some of said developed parameters 

to select a best frequency; and 
repeating the above steps with a new segment of data. 

15 

42. The method of claim 41, wherein said analyzing said frequency domain 
data to identify candidate peaks includes identifying all n candidate frequencies 
associated with the n largest amplitude peaks in power spectrum which are not harmonics 
of each other. 

20 

43. The method of claim 42, wherein n equals two. 



44. A circuit card for use in a pulse oximetry system to remove motion- 
induced noise artifacts from attenuated light signals, said circuit card comprising: 
25 a circuit board for mounting electronic circuitry and interfacing with the pulse oximetry 
system; 

a processor mounted on said circuit board for processing at least one input signal 

according to instructions; and 
memory for storing a computer program, wherein said memory is operably coupled to 
30 said processor, and wherein said computer program includes instructions for 
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implementing a method of removing motion artifacts from said attenuated light 
signals, said method comprising: 

acquiring a segment of red data and a segment of IR data from said attenuated 
light signals; 

5 conditioning said red segment of data and said IR segment of data for signal 

processing; 

transforming said conditioned segments of data into frequency domain data; 
analyzing said frequency domain data to identify candidate peaks; 
developing parameters associated with each of said candidate peaks; 
10 arbitrating between said candidate peaks using at least some of said developed 

parameters to select a best frequency; and 
repeating the above steps with a new segment of data. 

45. The circuit card of claim 44, wherein said processor is a digital signal 
15 processor. 

46. A pulse oximetry system for removing motion-induced noise artifacts 
from electrical signals representative of attenuated light signals comprising an input 
device, an output device, and a motion artifact circuit card, wherein said motion artifact 

20 circuit card comprises: 

a circuit board for mounting electronic circuitry and interfacing with the pulse oximetry 
system; 

a processor mounted on said circuit board for processing at least one input signal 
according to instructions; and 
25 memory for storing a computer program, wherein said memory is operably coupled to 
said processor, and wherein said computer program includes instructions for 
implementing a method of removing motion artifacts from said measured pulse 
oximetry signals, said method comprising: 

acquiring a segment of red data and a segment of IR data from each of the 
30 electrical signals representative of attenuated light signals; 
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transforming both said segment of red data and said segment of IR data into red 

and IR frequency domain data, respectively; 
identifying a plurality of candidate peaks from said red and IR frequency domain 

data; 

5 developing parameters associated with each of said plurality of candidate peaks; 

analyzing each of said plurality of candidate peaks with respect to said developed 
parameters; 

arbitrating between at least some of said plurality of candidate peaks using at least 
some of said developed parameters to select a best frequency; 
1 0 outputting pulse rate and saturation relating to said best frequency; and 

repeating the above steps for new segments of data. 

47. The pulse oximetry system of claim 46, wherein said processor is a digital 
signal processor. 

15 

48. A pulse oximetry system for removing motion-induced noise artifacts 
from electrical signals representative of attenuated light signals comprising an input 
device, an output device, and a motion artifact circuitry, wherein said motion artifact 
circuitry includes: 

20 a processor for processing at least one input signal according to instructions; and 

memory operably coupled to said processor for storing a computer program, wherein said 
computer program includes instructions for implementing a method of removing 
motion artifacts from said measured pulse oximetry signal, wherein said method 
comprises: 

25 acquiring a segment of red data and a segment of IR data from each of the 

electrical signals representative of attenuated light signals; 
transforming both said segment of red data and said segment of IR data into red 

and IR frequency domain data, respectively; 
identifying a plurality of candidate peaks from said red and IR frequency domain 
30 data; 
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developing parameters associated with each of said plurality of candidate peaks; 
analyzing each of said plurality of candidate peaks with respect to said developed 
parameters; 

arbitrating between at least some of said plurality of candidate peaks using at least 

some of said developed parameters to select a best frequency; 
outputting pulse rate and saturation relating to said best frequency; and 
repeating the above steps for new segments of data. 
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ABSTRACT OF THE DISCLOSURE 
A method for removing motion artifacts from devices for sensing bodily 
parameters and apparatus and system for effecting same. The method includes analyzing 
segments of measured data representing bodily parameters and possibly noise from 
5 motion artifacts. Each segment of measured data may correspond to a single light signal 
transmitted and detected after transmission or reflection through bodily tissue. Each data 
segment is frequency analyzed to determine up to three candidate peaks for further 
analysis. Each of the up to three candidate frequencies may be filtered and various 
parameters associated with each of the up to three candidate frequencies are calculated. 

10 The best frequency, if one exists, is determined by arbitrating the candidate frequencies 
using the calculated parameters according to predefined criteria. If a best frequency is 
found, a pulse rate and Sp0 2 may be output. If a best frequency is not found, other, 
conventional techniques for calculating pulse rate and Sp0 2 may be used. The above 
method may be applied to red and infrared pulse oximetry signals prior to calculating 

15 pulse rate and/or pulsatile blood oxygen concentration. Apparatus and systems disclosed 
are configured to perform methods disclosed according to the invention. 



20 
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